Density-functional theory study of half- metallic heterostructures: interstitial Mn in Si 
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Using density-functional theory within the generalized gradient approximation, we show that Si- 
based heterostructures with 1/4 layer (5-doping of interstitial Mn (Mni„i) are half- metallic. For Mni„t 
concentrations of 1/2 or 1 layer, the states induced in the band gap of 5-doped heterostructures still 
display high spin polarization, about 85% and 60%, respectively. The proposed heterostructures are 
more stable than previously assumed ^-layers of substitutional Mn. Contrary to wide-spread belief, 
the present study demonstrates that interstitial Mn can be utilized to tune the magnetic properties 
of Si, and thus provides a new clue for Si-based spintronics materials. 

PACS numbers: 75.50.Pp, 71.20.At, 75.70-i 
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Doping of semiconductors with magnetic 3d transi- 
tion metal impurities has attracted much interest: re- 
cent advances in spintronics have been largely due to 
magnetic semiconductors produced in this way. Their 
spin-dependent electric transport properties [ij could be 
exploited in novel magneto-electronic devices, e.g. spin 
filters. Since it has long been known that manganese im- 
purities can have a high magnetic moment [3,[3|, research 
has focused Mn-doped magnetic semiconductors based 
on GaAs and Ge 0, H, [^. Moreover, doping of GaAs 
by a (5-layer of Mn has been proposed theoretically 0], 
and an experimental realization has been reported re- 
cently Q. The calculations of spin-polarized transport 7] 
employing density-functional theory (DFT) within the 
local-spin-density approximation (LSDA) and pseudopo- 
tentials predicted that GaAs doped with a (5-layer of sub- 
stitutional Mn (Musub) on Ga sites is a two-dimensional 
(2D) ferromagnetic half-metal with large exchange cou- 
pling. Indeed, a Curie temperature of up to 250 K has 
been measured on Mn (5-doped GaAs samples Q . Adding 
magnetic functionality to the most common semiconduc- 
tor. Si, is still in its infancy, despite recent reports about 
ferromagnetism in Mn-doped Si created by ion implanta- 
tion Recently, Qian et al. [lo'| proposed, on the ba- 
sis of pseudopotential calculations within the generalized 
gradient approximation (GGA), that a Si heterostruc- 
ture with Mugub (^-doping should be a 2D ferromagnetic 
half-metal. So far, research on Mn-doped Si has concen- 
trated on substitutional Mn 11 , l2| , probably motivated 



by the physics of dilute magnetic semiconductors, where 
it has been argued that Mugub is crucial for ferromag- 
netism [1, [^, 14 1 . However, these substitutional Mn 



impurities in Si are energetically less stable than intersti 
tial ones [3, 3 iZ|- This has been considered a serious 
obstacle in creating magnetic Mn:Si calling for growth 
techniques far from equilibrium [l^. 

In this Letter, we investigate the role of interstitial Mn 
(Mnint) impurities for ferromagnetism in Si by means of 
all-electron full-potential DFT calculations. We propose 
a novel type of heterostructures with Muint (5-doping. 



For occupation of each forth interstitial site in a single 
layer by Mn, we find 2D half-metallic behavior with a 
Kohn-Sham gap of 0.5 eV in the electronic density of 
states in the majority spin channel. Even for higher dop- 
ing densities, the spin-polarization is still large, >60%, 
judged from the density of states at the Fermi level. 
This is different from the previously proposed Mugub 5- 
layers [lO|, where a (smaller) gap in the minority spin 
channel between Mn states of e and i2-type symmetry 
was responsible for the half-metallicity. Moreover, our 
proposed Mnint (5-layer is calculated to be 0.5eVperMn 
atom lower in energy than the heterostructure proposed 
by Qian et al. [lo| . in line with the known higher stabil- 
ity of isolated Mnjnt impurities compared to Mugub im- 
purities in Si 

0, 

Il7l |. Consequently, Si self-interstitials 
can destroy the half-metallicity of the heterostructure as- 
sumed by Qian et al. by energetically favorable site ex- 
change with Musub- 

Our proposed Si-based heterostructures are modeled 
by a Si(OOl) supercell with 16 atomic layers (see Fig. 1). 
Such thickness ensures that the central Si layer behaves 
as in the bulk, and the interaction between Mn layers is 
neghgible. Mn deposited on a Si(OOl) surface occupies 
preferably subsurface interstitial sites [l8[, and a c(4 x 2) 
arrangement for 1/4 monolayer (ML) Muint has been 
identified recently by a combined scanning tunneling mi- 
croscopy (STM) and DFT study [li|. Therefore we first 
study the c(4 x 2) heterostructure consisting of 16-layer 
Si and 1/4 ML Mnint in c{A: x 2) arrangement, as seen 
in Figs. 1(a) and 1(b). We also perform calculations for 
1/4 ML Mnint in p(2 x 2) arrangement (not shown). 
Secondly, we deal with a p{2 x 2) heterostructure with 
1/2 ML Mnint, as shown in Fig. 1(c). Finally we compare 
the (1x1) heterostructure having 1 ML Mni„t (Figs. 1(d) 
and 1(e)) with the (1x1) heterostructure having 1 ML 
Mnsub (Figs. 1(f) and 1(g)) studied by Qian et al. [l^]. 

We perform DFT-GGA calculations using the 
full-potential augmented plane-wave plus local-orbital 
method |2l| . which was shown to be appropriate for 
description of the Mn:Si system The lattice con- 
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FIG. 1: (a) (001)-2:y-plane view and (b) (110)-a;2:-plane view 
of the Si-based heterostructure (modeled by a 16-layer super- 
cell) with 1/4 layer of interstitial Mn in c(4 x 2) arrangement; 
(c) (OOl)-plane view of the heterostructure with 1/2 layer of 
interstitial Mn in p(2 x 2) arrangement; (d) (OOl)-plane view 
and (e) (llO)-plane view of the heterostructure with a full 
layer of interstitial Mn in (1 x 1) arrangement; and (f) (001)- 
plane view and (g) (llO)-plane view of the heterostructure 
with a full layer of substitutional Mn in (1 x 1) arrangement. 
Black balls show Mn and gray balls Si atoms. The silicon 
atoms in different layers are marked with indices. 

stant of 5.48 A for bulk Si calculated in GGA is used 
to construct the heterostructures we study. The muffin- 
tin sphere radii are chosen to be 1.11 A for both the 
Mn and Si atoms. Converged results are obtained at a 
cut-off energy of 13.8 Ryd for the interstitial plane-wave 
expansion, and with a set of 12x12x1 special k-points 
for integrations over the Brillouin zone of the (1x1) het- 
erostructure (and an equivalent density of k-points for 
larger supercells, e.g., 6x5x1 for the c(4 x 2) cell). All 
Si and Mn atoms except for the central Si layer are re- 
laxed until the calculated atomic force for each of them 
is smaller than 0.05 eV/A. 

The c(4 X 2) heterostructure (see Figs. 1(a) and 1(b)) 
displays a half-metallic electronic band structure as 
shown by the density of states (DOS) derived from the 
Kohn-Sham energy levels (see Fig. 2(a)). The spin-up 
channel (dashed blue lines) is insulating with a Kohn- 
Sham gap of 0.5 eV. The spin-down channel (solid red 
lines) is metallic, and the Fermi level lies 0.2 eV be- 
low the bottom of the spin-up conduction bands. While 
the Mn 3d spin-up orbitals are fully occupied, the pla- 



nar spin-down x'^ — orbitals form a conduction band 
with bandwidth of 0.7 eV, contributing to the DOS in 
the metallic channel. The narrower xz- and y2:-derived 
conduction bands are also partially occupied and add to 
the DOS at the Fermi energy The coplanar Si (Sio), 
the first-neighbor-layer Si (Sii) and the second-neighbor- 
layer Si (Si2) (cf. Fig. 1(b)) also contribute consider- 
ably to the half-metallic DOS through Mn 3d-Si 3s3p 
hybridizations. Such contribution decreases strongly for 
the third-neighbor-layer Si (Sia) and eventually vanishes 
for the central-layer Si (Sig). The Mn atom has a lo- 
cal spin magnetic moment (within the muffin-tin sphere) 
of 2.56 fj,B, and its four Si neighbors Sii (to which it is 
tetrahedrally coordinated) get spin-polarized with a spin 
moment of about 0.024 hb each. Any other Si atom has 
a spin moment less than 0.01 fiB- Taking also into ac- 
count the spin moment of 0.25 fj,B outside the muffin-tin 
spheres, stemming from the Mn 3d-Si 3s3p hybridized 
bands, the total spin moment is an integer of 3 ^b per 
supercell [l^]. Note that due to the long distance of 8.7 A 
between the center and corner Mn atom(s) in Fig. 1(a) 
the magnetic coupling turns out to be weak, but ferro- 
magnetic (FM): Our calculations show that the FM state 
is slightly more stable than the antiferromagnetic (AF) 
one by 5 meV/Mn, consistent with the short range of FM 
interactions in zincblende MnSi and MnGe [3, ■ 

The 1/4 ML p{2 x 2) heterostructure (not shown, but 
refer to Fig. 1(c) and remove the central Mn in the p(2 x 
2) cell) also displays half-metallicity, similar to the above 
discussed 1/4 ML c(4 x 2) heterostructure. The spin- 
down x^ — orbitals hybridized with Si sp and the spin- 
down xz and yz orbitals determine the half-metallic DOS 
(not shown). The p{2 x 2) superstructure is only sightly 
less stable than the above c(4 x 2) one by 20 meV per 
supercell, which shows that the c(4 x 2) superstructure 
gains a little more by lattice relaxation. 

For the 1 /2 ML p{2 x 2) heterostructure (see Fig. 1(c)), 
the DOS is displayed in Fig. 2(b), indicating a high spin- 
polarization at the Fermi level of 85%. The Mn spin 
moment is 2.72 ^b and the total one is 6.34 pB per su- 
percell (3.17 per Mn). Compared with the c(4 x 2) 
supercell, the conduction bands become broader due to 
decreased Mn-Mn distance and enhanced Mn d-Si sp hy- 
bridization. The FM exchange coupling of Mn to its four 
first neighbors at 5.5 A is quite large, J = 13 meV. This 
estimate is based on a Heisenberg model with local Mn 
spin of S = 3/2, and the calculated stability of the FM 
state over the AF state 475"^ = 120 meV per Mn atom. 
Thus the 1/2 ML p{2 x 2) heterostructure combines the 
advantage of a stronger in-plane FM coupling with a still 
high spin-polarization. 

Now we turn to the 1 ML (1x1) heterostructure (see 
Figs. 1(d) and 1(e)). Even in this case, spin-polarization 
is still > 60%, as shown by our calculations (Fig. 2(c)). 
The Mn and Si atoms form broad conduction bands. 
The Mn spin moment is 2.68 ^b and the total one is 
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FIG. 2: (Color online) DOS of (a) the 1/4-layer interstitial 
Mn c(4 X 2), (b) the 1/2-layer Mn p{2 x 2), and (c) the 1-layer 
Mn (1 X 1) heterostructures [cf. Figs. l(a)-l(e)]. The solid red 
(dashed blue) line refers to the down (up) spin component. 
The Fermi level is set to zero, (a) shows half-metallicity, (b) 
and (c) a high spin-polarization at the Fermi level. 



3.12 per (1 x 1) supercell. According to our calcu- 
lations, the Mn-Mn in-plane FM coupling turns out to 
be further increased: the FM state is more stable than 
the AF state by 570 meV/Mn. Since the direct mag- 
netic coupling between two neighboring Mn atoms (with 
distance of 3.9 A) is AF, we attribute the increased FM 
coupling with shorter Mn-Mn distance to an enhanced 
interaction of double-exchange type. It operates in the 
impurity band formed by the spin-down id states of Mn 
with t2 character, as seen in Fig. 2, and is mediated 
by spin-down Si sp states. As a result, we find the Si 
atoms to have small spin moments parallel to Mn in all 
the above calculations. 

Next we show that the heterostructure with 
Musub (Figs. 1(f) and 1(g)) proposed by Qian et 
al. [l^ is energetically unfavorable. Fig. 3(a) shows 
that our calculations yield half-metallic behavior with a 
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FIG. 3: (Color online) DOS of the 1-layer substitutional Mn. 
(a): (Ix 1) heterostructure shows half-metallicity; (b): p(2x2) 
heterostructure with extra 1/4-layer co- planar interstitial Si 
shows, after an energetically favorable Si-Mn site- interchange, 
the loss of half-metallicity. The solid red (dashed blue) line 
refers to the down (up) spin component. 



Kohn-Sham gap of 0.25 eV for 1 ML Mus^b (1 x 1), in 
agreement with Qian et al. [10]. The Mn spin moment 
is 2.97 fiB, each first-neighbor Si has a spin moment of 
— 0.116/is, and the total one is again integer, 3 /is per 
(1 x 1) supercell. However, the Mugub heterostructure 
turns out to be less stable than the Mujnt one by 
490 meV per (1 x 1) supercell, assuming bulk Si as 
reservoir for the extra Si. Knowledge about Mn-Si 
compounds corroborates the metastability of tetrahe- 
drally coordinated Mugub: Mn is known to prefer high 
coordination to Si, as exemplified by the seven-fold 
Mn-Si coordination in the ground-state crystal structure 
of MnSi, or the eight-fold coordination in the CsCl 
structure [ll|. In contrast, the hypothetic zinc-blende 
structure of MnSi with tetrahedral coordination of Mn to 
Si is calculated to be extremely unstable, 2.3 eV higher 
in energy than the ground state. The heterostructures 
with Muint proposed by us, having four second-neighbor 
Si atoms (2.7 A) in addition to four first-neighbor Si 
atoms (2.4 A), comply with the rules for the optimal 
Mn-Si coordination derivable from bulk phase stability. 

We note that both types of (5-layers considered are 
only metastable in the thermodynamic sense, i.e. they 
could be destroyed by annealing. However, the ten- 
dency of Mn to form ordered subsurface layers, as ob- 
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FIG. 4: Schematic view of the electronic states of the (a) 
interstitial and (b) substitutional Mn in Si (2^. The filled, 
patterned, and open circles indicate full, partial, and zero 
occupation, respectively. Both the size of the circles and the 
numbers next to the levels indicate their occupation number. 
The half-metallic DOS results in (a) from the band, but 
in (b) from the t2i antibonding band. 



served recently [19|, suggests that growth kinetics could 
be exploited to fabricate the Mnint heterostructure pro- 
posed by us. Although the reported Mn-induced c(4 x 2) 
structure has so far been observed only in islands, ex- 
tended layers might be obtainable by improved growth 
techniques, and could be overgrown at lower tempera- 
tures with a capping layer. After successful preparation, 
the heterostructures are still jeopardized by degradation 
due to highly mobile [2^ Si self-interstitials: By study- 
ing the 1 ML Mnsub p(2 x 2) heterostructure with an 
extra 1/4 ML co-planar interstitial Si, we find that in- 
terstitial Si is capable of site exchange with one Musub, 
thereby considerably lowering its energy by 1.2 eV. This 
irreversibly destroys the half-metallicity of the Mugub- 
heterostructure (cf. Ref. (lOj) by degrading the spin- 
polarization at the Fermi level to about 40%, as seen 
in Fig. 3(b). In the Muint -heterostructure, trapped co- 
planar Si intcrstitials have a less drastic effect on the 
electronic structure, but they could shift the Fermi level 
out of the half-metallic gap if present in high concentra- 
tions. 

We stress that the half-metallicity in our proposed 
Mnint-beterostructure (metallic minority spin channel) 
and in the Musub-heterostructure of Qian et al. (metal- 
lic majority spin channel) have different physical origin. 
This is explained in the level diagram of Fig. 4. The 
Muint 3d states are energetically located near the top of 
the Si-host valence bands (Fig. 2(a)). The weak crystal 
field with Td symmetry and strong Hund exchange stabi- 
lize a high-spin state for Mnint [SS]- The 4s -^3d elec- 
tron transfer leads to a 3ci'' [t^^e^tl^) configuration (Fig. 
4(a)). Thus the partially filled t2i band is responsible for 
the metallic spin-down DOS seen in Fig. 2(a). For Mugub, 
however, the Mn states of t2 symmetry mix strongly with 
Si sp states of the same symmetry, thus splitting the t2 
states into a bonding and an antibonding band, lying in 
the valence and in the conduction band of the Si host. 



respectively (cf. Fig. 3(a)). The Mugub 3d ^2"^^ ^'^d 
As^y , in total four electrons, fill up the bonding states, to- 
gether with the Si sp electrons (Fig. 4(b)). Note that the 
4s — > 3d electron transfer in the Mngub is not as complete 
as in the Muint, as our calculations find the Mngub 3d 
occupation (within the mufRn-tin sphere) to be less than 
the Muint one by 0.2 electrons. Since the non-bonding 
e| states (lying below the t2\ states for Td symmetry) are 
fully occupied by two electrons, the remaining one elec- 
tron resides in the t2] antibonding band that gives rise 
to the DOS in the spin-up channel in Fig. 3(a). 

To conclude, we proposed a stable Si-based het- 
erostructure using (5-doping by interstitial Mn, and 
proved it to be half-metallic for 1/4 monolayer of Mn, 
using DFT calculations. A recent combined STM-DFT 
study [Ig] lets us expect that such structures could be 
prepared by molecular beam epitaxy. In order to achieve 
ferromagnetic ordering at room temperature, one should 
aim at higher Mn concentrations. Even for a full (5-layer 
of Muint , the spin polarization of the conduction electrons 
should still be as high as 60%. 
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